Astrophysical Measurement of the Equation of State of Neutron Star Matter 



Feryal Ozel, 1 Gordon Baym, 2 and Tolga Giiver 1 

1 University of Arizona, Department of Astronomy and Steward Observatory, 933 N. Cherry Ave., Tucson, AZ 85721 
2 Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL 61801 

We present the first astrophysical measurement of the pressure of cold matter above nuclear 
saturation density, based on recently determined masses and radii of three neutron stars. The 
pressure at higher densities are below the predictions of equations of state that account only for 
nucleonic degrees of freedom, and thus present a challenge to the microscopic theory of neutron star 
matter. 
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Neutron stars probe the dense QCD phase diagram at 
lower temperatures and higher baryon densities, in con- 
trast to the higher temperature-lower density regime in 
the earlyuniverse and in ultrarelativistic heavy-ion colli- 
sions [J1H • The baryon density, p, in stellar interiors may 
reach an order of magnitude beyond nuclear saturation 
density, p ns ~ 2.7 x 10 14 g cm~ 3 ~ 0.16 fm~ 3 , and cross 
into a regime where quark degrees of freedom are excited 
or matter undergoes a meson condensation phase transi- 
tion. The pressure of matter at these high densities can, 
as we show here, be extracted from current neutron star 
mass-radius determinations Q, and crucially constrains 
calculations of high density neutron star matter. 

The equation of state (EoS) of supranuclear matter 
determines the dividing line between neutron stars and 
black holes, and directly impacts the mechanism as well 
as outcomes of supernova explosions and the numbers of 
neutron stars and black holes in the Galaxy [1, [||. In 
certain models, it affects the mechanism and duration of 
gamma-ray bursts [6]. Accurate evolutions of inspiral- 
ing neutron-star binaries and the collapse to black holes, 
needed to calculate gravitational wave signals, depend 
sensitively on the assumed EoS 

Microscopic calculations of the EoS of neutron-star 
matter have been based on a variety of inputs. The ap- 
proach most firmly founded on experiment in the region 
of p ns is to determine two-body potentials from nucleon- 
nucleon scattering data below 350 MeV and properties 
of light nuclei, supplemented by a three-body poten- 
tial [1, Hj]. Such calculations, accurate in the neighbor- 
hood of p nsi have fundamental limitations. Beyond a few 
times p ns the forces between particles can no longer be 
described via static few-body potentials; since the char- 
acteristic range of the nuclear forces is ~ l/2m„, where 
m„ is the mass of the pion, the parameter measuring the 
relative importance of three and higher body forces is 
~ p/(2m 7r ) 3 ~ 0.35p/p ns . Thus, at p ;§> p ns a well defined 
expansion in terms of two-, three-, or more, body forces 
no longer exists. EoS based on nucleons alone do not take 
into account the rich variety of hadronic degrees of free- 
dom that enter with increasing density. In addition, pion 
condensates [1, [l(| [H| or kaon condensates [HI, can 
enter at higher densities. Field-theoretic models based 
on nucleons interacting via meson exchange include, e.g., 
Ref [HJ|; see Ref. 15 1 for a general summary of EoS. How- 



ever, one cannot assume that matter at higher densities 
can even be described in terms of well-defined "asymp- 
totic" laboratory particles. More realistically, one ex- 
pects in dense matter a gradual onset of quark degrees 
of freedom, not accounted for by nucleons interacting via 
static potentials. Indeed once nucleons overlap consider- 
ably, the matter should percolate, opening the possibility 
of their quark constituents propagating throughout the 
system. Such additional degrees of freedom should lead 
to softening of the EoS, consistent with our findings here, 
and thus to a lower maximum neutron star mass. Ow- 
ing to difficulties of determining the neutron star matter 
EoS from first principles, neutron star observations be- 
come an important input in determining the EoS at high 
density, and constraining microscopic calculations. 

A number of astrophysical observations have focused 
on measuring neutron star radii, R, and masses, M, in 
an attempt to constrain the uncertainties in the EoS 
(summarized in Ref. [Hj]). The measurement of post- 
Newtonian parameters of double neutron stars provide 
precise determinations of their masses [l||. Glitches ob- 
served in radio pulsars lead to constraints on the mo- 
ment of inertia, and therefore, on neutron star masses 
and radii [ItJ . Observations of the thermal emission from 
accreting neutron stars in quiescence and from millisec- 
ond X-ray pulsars result in broad, correlated constraints 
on neutron star masses and radii [l8| . 

Recent advances in both astrophysical techniques and 
neutron-star atmosphere modeling allow us for the first 
time to measure the pressure of neutron star matter at 
supranuclear densities directly from observations. Refer- 
ences [15j suggested that the radii of neutron stars are 
a good indicator of the pressure at roughly twice the 
nuclear saturation density. This argument was further 
extended in [ijj [2(| , which demonstrated that knowing 
the properties of ultradense matter at three fiducial den- 
sities allows one to reproduce macroscopic neutron star 
properties, including the mass-radius relation and stellar 
moment of inertia. Conversely, three distinct measure- 
ments of neutron star masses and radii, as we use here, 
is sufficient to infer a piecewise polytropic EoS of matter 
at supranuclear densities [19j. 

Observations of multiple spectroscopic phenomena 
during thermonuclear bursts from X-ray binaries have 
enabled the tightest measurements neutron star radii 
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and masses to date [3J- The long-term monitoring of 
burst sources with the Rossi X-ray Timing Explorer, with 
its excellent photon statistics, has resulted in a large 
(> 1000) database of bursts [2l|, from which systematic 
uncertainties can be determined and controlled, and ideal 
sources that act as standard candles can be identified. 
High resolution X-ray spectroscopy with the Chandra X- 
ray Observatory and XMM-Newton has led to a detailed 
measurement of the soft X-ray spectra of bursters and re- 
duced the uncertainties introduced by interstellar extinc- 
tion [22j]. Finally, pointed optical/infrared observations 
with the Hubble Space Telescope and large ground-based 
facilities (e.g., the Magellan telescope) have substantially 
improved distance measurements to these sources (see 
Ref. Q and references therein). 

The wealth of such high quality data allows us to em- 
ploy a novel approach, combining different spectroscopic 
measurements to break the degeneracies between neutron 
star masses and radii inherent to each observable |23u24j. 
The first observable is the apparent surface area during 
the cooling phase of the bursts, 

R 2 / 2GMY 1 

where D is the distance to the source, and f c is a cal- 
culated ratio between the spectral (color) and effective 
temperature, T, of the emerging radiation that accounts 
for the non-Planckian spectrum of the burst. Because 
the emitted luminosity is cx T 4 , the apparent surface area 
in Eq. ([T]) shows the same T dependence, which we ab- 
sorb into the definition of f c . The apparent surface area 
remains constant in time and is highly reproducible in 
multiple events from the same source, indicating that the 
entire neutron star surface, rather than a variable area 
on the surface, participates in the burst emission. 

The second phenomenon occurs in a subset of bursts, 
when the flux becomes so high that it exceeds the local 
Eddington limit and lifts the photosphere of the neutron 
star. The flux achieved during these events is also highly 
reproducible for a large number of sources including the 
three discussed below, and is related to the neutron star 
mass and radius through 

_GMc/ 2GAA 1/2 

This "touchdown" flux is evaluated at the moment when 
the photosphere has receded back to the neutron star 
surface (see also Ref. Q). In the above equation, fc os = 
0.2(1 + X) cm 2 g -1 is the electron scattering opacity in 
the stellar atmosphere. In the atmospheric models, we 
considered a wide range of hydrogen mass fraction X 
consistent with the properties of each binary system. 

Combining the measurements of A and F with the dis- 
tance D to each source, we obtain tight, uncorrelated 
constraints on the masses and radii of neutron stars. We 
have applied this technique to three sources, the neu- 
tron stars in the binaries 4U 1608-248, EXO 1745-248, 
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FIG. 1: The 1- and 2-<r confidence contours for the masses 
and radii of three neutron stars in the binaries 4U 1608—248 
(green/red), EXO 1745-248 (yellow/blue), and 4U 1820-30 
(cyan/magenta), compared with predictions of representative 
EoS (see text for details). The details of the measurements are 
described in Refs. 0|. The diagonal lines are the black-hole 
event horizon (solid) and Buchdahl (dashed) 25] limits. 



and 4U 1820-30, and show in Fig. [T] the 1- and 2-cr 
confidence contours of their masses and radii determi- 
nations 3]. The results are a set of uncorrelated mea- 
surements of neutron star masses and radii. (An earlier 
measurement of the mass and radius of the neutron star 
in EXO 0748—676 [23| was based on the identification of 
atomic line features in its X-ray spectrum with gravita- 
tionally redshifted lines from its surface, which has since 
been shown to be inconsistent (27j with the recent mea- 
surement of its rapid spin frequency [28j]). 

The measurements in Fig. [TJ incorporate the correc- 
tions and systematic uncertainties associated with the 
modeling of emission from the hot surfaces of neutron 
stars following a thermonuclear burst, the composition 
of the neutron star surface, as well as statistical or sys- 
tematic uncertainties in the distances to the binaries. In 
all three measurements, uncertainties arising from sub- 
traction of the background flux are negligible, because in 
each source, the luminosity from the neutron star surface 
exceeds the accretion luminosity by more than a factor 
of ten. Since not every combination of observables leads 
to a solution for M and R, we converted the probability 
densities over the measured fluxes, apparent areas, and 
distances to those over the neutron star mass and radius 
following standard Bayesian statistics Q. As a result, the 
uncertainties in the mass and radius, of order 15%, are 
smaller than those of the individual spectroscopic quan- 
tities. Note that the confidence contours shown in Fig. Q] 
correct a small numerical error in the Jacobian transfor- 
mation of Refs. Q. 

As noted, three distinct measurements of neutron star 



3 



E 

4— 

> 



o 



> 



o 




log(P 2 /MeV fm ) 

FIG. 2: The pressure of cold matter at (top) 7.4 and 3.7 p„ 
and (bottom) 1.85 and 3.7 p nB . 



masses and radii allow us to infer a piecewise EoS of mat- 
ter at supranuclear densities This approach makes 
the explicit assumption that the density in the neutron 
star surface layers smoothly reduces to zero, and is, there- 
fore, not applicable to strange quark matter stars that are 
not gravitationally bound. In particular, we follow the 
procedure in Ref. [19( to convert these measurements to 
probability densities over the pressures of neutron star 
matter at three fiducial baryon densities p\ = 1.85p ns , 
p% = 3.7pns, and p^ = 7.4p ns . In these calculations, we 
include the full probability density for all three sources. 
We supplement this procedure with the requirement of 
causality, rejecting combinations of pressures for which 
the sound speed is larger than the speed of light [l9j . 
Fig. [U shows the confidence contours of different pairs of 
pressures (integrated over the third pressure). 

The pressure at p\ is only weakly constrained because 
such low densities are important only in determining 
the macroscopic properties of neutron stars with masses 
smaller than those in our measured sample [l9j . In con- 
trast, the pressures at pi and ps are constrained to within 
a factor of - 10 3 . 

In detail, we describe these data with a phenomeno- 
logical piecewise polytropic EoS above density po — 
10 14 g cm -3 , which we smoothly connect to the SLy 
EoS at lower densities [26]. Between p and p4 we fit 
the energy density in the interval pi-i < p < pi as 



TABLE I: 

Measured pressures at three supranuclear densities, in 
MeV fm -3 , together with Po taken from the low density SLy 
EoS [2f|. 
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FIG. 3: Pressure vs. baryon density at the three points, Pi, 
P%, and P3, together with the fitted EoS The shaded 

region shows the uncertainties in the determination. 



e = otip + /3ip Ti , from which we derive the pressure, 
P = p 2 d(t/p)/dp = (Ti - l)Ap r *. Since the Tolman- 
Oppenhcimer-Volkoff equation relates e and P to the 
mass and radius of the star, one only can determine the 
baryon density p from e and P data to within a scale fac- 
tor; to determine the scale we connect our fit to the tab- 
ulated low density SLy EoS. For < p < Pi, the effec- 
tive polytropic index is r. i = log(Pi/i\-_i)/log(/9i//3i_i), 
the pressure is 



P = Pi 
and in the energy density 



Pi-i (Ti - l)pi-i 



Pi-i 
Pi 



(3) 



(4) 



Table 1 shows the fitting parameters. For p > p%, we 
extrapolate the last polytropic relation. 

The data present a clear challenge to microscopic nu- 
clear calculations. Figs. [5] and [3] compare the best fit 
values of the pressures at the three fiducial densities 
with those predicted by a representative sample of EoS 
based on a wide range of input physics and computational 
methods, from nucleonic: variational chain summation 
with the AV18 potential, UIX three-body potential plus 
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relativistic boost corrections AP4 0; Dirac-Brueckner- 
Hartree-Fock MPA1 ; and relativistic mean fields MSI 
(l4j . plus kaons GS1 [13[. (Although our procedure is not 
applicable to the u,d,s quark matter SQM1 [la ], its pre- 
dictions can nevertheless be compared directly to the in- 
dividual mass and radius measurements shown in Fig. [I}. 
Our measurements clearly discriminate between different 
predictions, and indicate that the EoS based on nucleons 
alone, AP4, MP1 and MSI, are too stiff at higher density 
- a conclusion also borne out by the comparisons of M 
vs. R in Fig. [T] which show that the predicted radii at 
the measured masses are too large. As one sees clearly 
in Fig. [3l the data call for a softer EoS, as would be pro- 
duced by including degrees of freedom beyond nucleons, 
e.g., hyperons, mesons, and quarks, or possibly produced 
by a better description of nucleonic interactions. In par- 
ticular a softer EoS, by allowing larger central densities 
than an EoS with only two and three body nucleonic in- 
teractions (e.g., AP4 with a central density ~ 7p ns at the 
maximum mass ~ 2.2M@) would, in the naive picture of 
a sharp phase transition from nucleonic to quark matter, 
allow quark cores to appear in more massive stars. 
The allowed range of pressures seen in Fig. [3] also 



reduces the maximum neutron star mass compared to 
that predicted by purely nucleonic EoS. It is, never- 
theless, consistent with the recent measurement of a 
1.97 ± 0.04 Mq neutron star by observations of Shapiro 
delay [29|. 

Our measurements can be confirmed or tested by ad- 
ditional mass-radius measurements in bursting sources, 
improved distance measurements by space-based interfer- 
ometers, and more definitively by observations of other 
phenomena that probe the masses and radii of neutron 
stars such as gravitationally redshifted absorption lines 
or flux oscillations that depend on surface gravity. 
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